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The presence of dynamical heterogeneities, i.e. 
nanometer-scale regions containing molecules re- 
arranging cooperatively at very different rates 
compared to the bulk ]§], Q|, is increasingly be- 
ing recognized as crucial in our understanding of 
the glass transition, from the non-exponential na- 
ture of relaxation, to the divergence of the re- 
laxation times [5|. Recently, dynamical hetero- 
geneities have been directly observed experimen- 
tally @, @, H, [|| and in simulations [l3| . However 
a clear physical picture for the origin of these het- 
erogeneities is still lacking. Here we investigate 
a possible physical mechanism for the origin of 
dynamical heterogeneities in the non-equilibrium 
dynamics of structural glasses. We test the pre- 
dictions regarding universal scaling of fluctua- 
tions derived from this mechanism against sim- 
ulation results in a simple binary Lennard-Jones 
glass model, and find that to a first approxima- 
tion they are satisfied. We also propose to apply 
the same kind of analysis to experimental data 
from confocal microscopy in colloidal glasses. 

Supercooled liquids approaching the glass transition 
display increasingly slow dynamics, until eventually they 
cannot equilibrate in laboratory timescales [l[ . One con- 
sequence of this fact is physical aging, i.e. the break- 
down of time translation invariance (TTI): the correla- 
tion C(t,t w ) between spontaneous fluctuations of an ob- 
servable at times t (the final time) and t w (the wait- 
ing time) are nontrivial functions of t and t Wl as op- 
posed to being functions of the time difference t — t w . In 
many cases, the two-time correlation C(t,t w ) in an ag- 
ing system separates into a fast, time translation invari- 
ant contribution Cf as t(i — t w ), and a slow contribution 

Cslow 

(t,t w ) 0]: C(t,t w ) — Cf as t(i — t w ) + C s i ow (t,t w ). 
In the case of a structural glass, two-step relaxation is 
observed: the fast term corresponds to localized fluctu- 
ations of individual particles inside their cages, and the 
slow term corresponds to longer time scales, in which 
cages break down and the system structurally relaxes. 
For some systems, the slow part of the correlation has 
the form 0] C s i ow (t,t w ) = C s \ ow (h(t)/h{t w )), where h(t) 
is some monotonically increasing function. For example, 
in the case of domain growth, h(t) is proportional to the 
domain size Q. 

Recently, it has been proven that, in the limit of long 



times, the dynamics of a class of spin-glass models is 
invariant under global reparametrizations t — > h(t) of 
the time [lOj. This result has been used to predict 
the existence of a Goldstone mode in the nonequilib- 
rium dynamics, associated with smoothly varying lo- 
cal fluctuations in the reparametrization of the time 
t — -> h r (t) = e Vr ^ t \ These fluctuations have been phys- 
ically interpreted to represent local fluctuations of the 
age of the sample 
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12j . In the cases where the global 
two-time correlation exhibits h(t) / h(t w ) scaling, a simple 
Landau theory approximation for the dynamical action 



predicts [Hi[l2, 14 1 that the full probability distribution 
p(C r (t, t w )) of local correlations C r (t, t w ) depends only on 
the global correlation C g \ Q h a \{t,t w ). In all of this discus- 
sion, only the slow part of the correlation is considered, 
and any effects due to the fast part of the dynamics are 
neglected. 

In this Letter, we examine whether this theoretical 
picture can be extended to explain dynamical hetero- 
geneities in the aging of structural glasses, by present- 
ing the first detailed characterization of the behavior of 
fluctuations in the aging of a continuous-space, quasi- 
realistic structural glass model. The presence of a global 
symmetry under time reparametrization in the generat- 
ing functional for the dynamics has been proven only for 
a class of spin glass models, and it is far from obvious 
whether it applies in off-lattice, quasi-realistic models, 
describing structural glasses. Addtionally, the question 
about the aging behavior of local observables in struc- 
tural glasses remains mostly open. Almost all simulations 
in glass forming liquids have focused on the (equilibrium) 
supercooled liquid [ijj], although there are some recent 
results for kinetically constrained spin models of glassy- 
ness 14j, and an earlier work exploring dynamic spatial 



correlations for soft spheres in the aging regime [15]. 

We probe individual particle displacements along one 
direction Axj(t, t w ) — xj(t) — Xj(t w ) (where j is the par- 
ticle index) ; and also local, coarse-grained two-time func- 
tions: the correlator 



C r (t, t w ) 



1 



N(B r ) ^ 



cos(q-(r^)-r,(^))), (1) 



and the mean square displacement 
1 



A r (t,t w ) = 



N(B r ) 



E ( r j(*)- r i(U) 2 



(2) 



i(t w )<=B r 
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Here we consider a coarse graining cubic shaped box B T of 
side £ around the point r in the system, and the sums run 
over the N(B r ) particles present at the waiting time t w in 
the box B r . We choose a value of q that corresponds to 
the main peak in the structure factor S(q) of the system, 
q = 7.2 in Lennard- Jones units. 

These definitions are inspired by the analogous defini- 
tions in the case of spin glasses [ll], [l2[ , and can be ap- 
plied both to analyze data obtained from simulations and 
from confocal microscopy experiments. The global quan- 
tities C g iobai(i, tw) (incoherent part of the intermediate 
scattering function) and A g i bai(i) t w ) (mean square dis- 
placement) are defined by extending the sum to the whole 
system in Eq. (p} and Eq. ([2]) respectively. 

We performed 250 independent molecular dynam- 
ics runs for the binary Lennard- Jones (LJ) system of 
Ref. [loT ] , which has a mode coupling critical temperature 
T c = 0.435. A system of 8000 particles was equilibrated 
at a temperature To = 5.0, then instantly quenched to 
T = 0.4, and finally it was allowed to evolve for 10 5 
Lennard-Jones time units. The origin of times was taken 
at the instant of the quench. 

As mentioned above, we want to test whether the prob- 
ability distributions of local correlations depend on the 
two times t, t w only through the values of a global corre- 
lation. Unlike in spin systems, in the case of a struc- 
tural glass there are different global correlation func- 
tions associated with the microscopic degrees of freedom 
in the system, including the self-intermediate scattering 
function C K \ \ >a \{t,t w ) and the mean square displacement 
A g iobai(i, tw). This leads us to testing two possible hy- 
pothesis: a) that the probability distributions coincide 
when C g iobai(i, tw) is kept constant; and b) that the prob- 
ability distributions coincide when A g i b a i(t, t w ) is kept 
constant. Extensive checking of our numerical data indi- 
cates that a) is approximately satisfied by the probability 
distributions of C r (t, t w ), A T (t,t w ) and Ax( t, t w ), but b) 
is satisfied only to a much lesser degree [171 ]. 

In Fig. [T] we present our results for the probability 
distribution p(C r (t,t w )) of the local intermediate scat- 
tering function for waiting times t w — 30.20, • • • , 30200, 
and final times t chosen so that C g iobal(i) t w ) G 
{0.1,0.3,0.5,0.7}. We observe that the data approxi- 
mately collapse for each value of C g i bai(i, t w ). This kind 
of collapse is also observed in simulations in a 3D spin 
glass model, but in the case of the spin glass model, the 
collapse is more precise than here. Unlike the case of 
the 3D spin glass model, the position of the peak in the 
distribution p(C r ) is strongly dependent on the value of 
Cgiobai(i, tw)- The skewness of p(C r ) also depends dra- 
matically on C g iobai(£, t w ): it goes from highly skewed for 
Cgiobal(t) t w ) — 0.7 to almost symmetric around its peak 
for Cgiobai^, tw) = 0.1. As shown in the figure, the dis- 
tributions p(C r (t,t w j) for C g iobai(i, t w ) — 0.1 can be well 
approximated by a Gaussian fit, but this does not hold 
for larger values of C g iobal(£) t w ). Unlike previous results 
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FIG. 1: Probability distributions p(C r (t,t w )) for values t w 
in the range 30.20 - 30200 (as indicated by the key in the 
figure), plotted for final times t chosen so that C g i b a l(i, tw) € 
{0.1,0.3,0.5,0.7}. Coarse graining box size I w 0.11L, (with 
L = linear size of the simulation box). The curves collapse 
into four groups, corresponding to G g i b a l = 0.7,0.5,0.3,0.1 
(ordered from highest to lowest value of CV at the peak). A 
gaussian fit to the data for C = 0.1 is also shown. Top panel: 
linear scale. Bottom panel: logarithmic scale. 



in kinctically constrained models [14j , in the present case, 
Gumbel distributions do not provide good fits to p{C r ). 
To characterize the weak dependence of the probability 
distributions on waiting time at fixed C g i bai(i, t w ), in the 
top panel of Fig. [3] we plot the centered second moment 
of the distributions p(C r ) as a function of waiting time, 
for fixed C g iobal(M«;) S {0.1,0.3,0.5,0.7}. The depen- 
dence on t w is so weak that both a logarithmic form and 
a power law form (with powers in the range 0.01 — 0.07) 
provide a good fit. 

In Fig. [5] we present our results for the probability 
distribution p(Ax(t, t w )) of the particle displacements 
Axj(t, t w ) = x j{t) — Xj{tw) along one direction. In the 
top panel of Fig. [5] we can observe that these data also 
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FIG. 2: Top panel: Probability distributions p(Ax(t,t w )) 
for values t w in the range 30.20 - 30200 (as indicated by the 
key in the figure), plotted for final times t chosen so that 
Cgiobal(i) tw) £ {0.1,0.3,0.5,0.7}. The curves collapse into 
four groups, corresponding to C g i b a i = 0.1,0.3,0.5,0.7 (or- 
dered from lowest to highest value of p(Ax) at the peak). 
Bottom panel: Tails of the distributions p(Ax(t,t m )), for 
C = 0.5 and t w = 30.20,302,3020,30200 (from narrower 
to wider tail). Symbols: results from simulation. Lines: 
fits to the data for \Ax\ > 0.5 using a stretched exponen- 
tial form p(Ax(t,t w )) ~ Af exp( — \Ax/af ), with exponents 
(3 = 1.11, 1.01, 0.90, 0.81 respectively. 



approximately collapse for each value of C g \ bai(t,t w ). 
In the bottom panel we have a closer look at the tails 
of p(Ax(t,t w )). We find that the distribution is non- 
gaussian, as was observed in experiments in colloidal 
glasses in the supercooled regime 0- The tails of the 
distribution can be fitted with a stretched exponential 
form p(Ax) « A/"exp(— | Ax/af), and they become more 
prominent as t w grows (for constant C g i bai(i, t w )). In- 
deed, as shown in the bottom panel of Fig. [3J the expo- 
nent (3 decreases from values above unity at short waiting 
times t w to values of around 0.8 at much longer t w . 



FIG. 3: Evolution of the probability distributions, as a 
function of the waiting time t w , at constant C g i b a i(i, t w ) £ 
{0.1, 0.3, 0.5, 0.7}. Top panel: second moment of p(CV(t, i™)), 
together with fits to the functional forms: mo(t w ) a (full lines) 
and m \og(t m /to) (dotted lines). Bottom panel: stretching 
exponent /3 for the tails of p( Ax), as a function of the waiting 
time t m , at constant C g i b a i(i, tw) £ {0.1,0.3,0.5,0.7} (the 
lines are only guides to the eye). 



To summarize, in this Letter we have presented the 
first detailed characterization of non-equilibrium fluctu- 
ations in the aging regime in a continuous-space, quasi- 
realistic structural glass model. It is also the first test for 
the possible presence in such models of a Goldstone mode 
associated with reparamctrizations of the time variable, 
which could explain the behavior of fluctuations in the 
aging regime. We have found that, as expected from the 
Goldstone mode picture, the probability distributions for 
the local fluctuating two-time quantities are, to a first 
approximation, invariant when the global intermediate 
scattering function C g i bai(i, t w ) is kept constant. 

Besides the scaling predicted by the Goldstone mode 
picture, a detailed analysis of the probability distri- 
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butions for local observables uncovers that: i) unlike 
the case of spin glasses, the position of the peak for 
the probability distributions of the local correlation 
p(C r (t, t w )) shifts dramatically as a function of the value 
of Cgiobai (t, t w ) ; ii) the distribution p(C r (t,t w )) evolves 
gradually from being highly skewed and non-gaussian for 
larger values of C g i bai(i, t w ) to being unskewed and very 
close to gaussian for small values of C g i b a i(i, t w ); iii) 
a more detailed examination of the probability distri- 
butions at constant C g i t>ai(t,t w ) reveals that the mo- 
ments of the distributions evolve smoothly as a function 
of the waiting time t w , without displaying any obvious 
characteristic timescale (this may be due in part to the 
fact that the dynamic correlation length in the system 
is growing as a function of t w , but for the timescales of 
the simulation it is not yet larger than the size of the 
coarse graining box used [ijj], leading to a gradual in- 
crease in the variance of all coarse grained quantities) ; 
iv) the probability distributions of one-dimensional dis- 
placements p(Ax(t,t w )) are clearly non-gaussian, as in 
confocal microscopy experiments in the supercooled liq- 
uid regime 0, and their tails can be fit by stretched 
exponential forms, with exponents [3 that decrease from 
/3 > 1 to [3 w 0.8 as t w increases (for fixed C g i bai(i, t w )). 

We conclude that the probability distributions of local 
observables in the aging regime are consistent with the 
presence of a Goldstone mode controlling the nonequilib- 
rium dynamics of a structural glass. We suggest that a 
direct experimental test of this picture could be provided 
by applying a similar analysis to experimental data from 
confocal microscopy in colloidal glass systems. 
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